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Abstract 

We consider a (1+1) dimensional scalar field theory that supports oscillons, which are localized, 
oscillatory, stable solutions to nonlinear equations of motion. We study this theory in an expanding 
background and show that oscillons now lose energy, but at a rate that is exponentially small when 
the expansion rate is slow. We also show numerically that a universe that starts with (almost) 
thermal initial conditions will cool to a final state where a significant fraction of the energy of the 
universe — on the order of 50% — is stored in oscillons. If this phenomenon persists in realistic 
models, oscillons may have cosmological consequences. 

PACS numbers: 11.27.+d 05.70.Ln 98.80.Cq 
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INTRODUCTION 



A wide range of nonlinear field theories have been found to contain long lived, localized, 
oscillatory solutions to their equations of motion, known as oscillons or breathers. In the 
best-known examples [l|, [2|, conserved charges guarantee the existence of exact periodic 
solutions. However, in many cases where such arguments are not available, objects with 
similar properties have also been observed [sl. 13. Isl. IgTItI. H. [qI. [lO . 11, 12, 13, 14, 15, 16, 13, 18 . 

In this work we consider analytically the evolution of an oscillon in an inflating back- 
ground. We find that the oscillon is no longer localized and contains a "tail" that slowly 
leaks energy away. However this decay rate is exponentially suppressed, so we still expect 
there to exist long-lived objects. This result is in qualitative agreement with numerical work 
in [1^], which studied the long-term evolution of a single oscillon in an expanding universe 
background and found that oscillons remain stable for an exponentially long time provided 
that the horizon is far larger than the width of an oscillon. 

Once a model has been shown to contain oscillons, it is natural to ask how easy it is for 
these coherent objects to form from generic initial conditions. Ref. [13] showed that oscillons 
can emerge from a rapid "quench" in which the background potential is suddenly changed, 
throwing the system far out of equilibrium. Here we consider the opposite situation: We 
begin with an (almost) thermal distribution at high temperature and gradually cool the 
system by coupling it to an expanding background. This setup is suggestive of a situation 
that could arise in the early universe, for example as the universe cools after reheating or a 
phase transition. If stable oscillons formed in these situations, they mi ght have cosmological 



consequences: e.g., see [18[. Although it did not consider oscillons, [20[ found significant 
non-thermal effects in electroweak baryon number violating processes. Ref. [2l| studied 
a case similar to ours, but involving a different type of oscillon that is stable only if its 
amplitude exceeds a certain critical value. This situation is quite different from our model, 
in which there is no such threshold amplitude. 

For numerical convenience, we work with a single scalar field in one dimension, though 
we have seen qualitatively similar results in simulations of the two-dimensional scalar model 
of [gI in an expanding background. We start at a temperature for which the universe is 
dominated by radiation and allow the universe to expand until it contains only a cold, 
pressureless dust of both oscillons and fundamental excitations of the field. A late-time 
snapshot of the energy density of such a configuration as a function of position is shown in 
Figure [H where the sharp spikes are oscillons. We find that a sizeable fraction of the energy 
of this final state — of order 50% or more — is stored in oscillons. This result persists even 
for very small coupling constants, where quantum effects are small and do not affect our 
classical field theory analysis. 



MODEL 

We consider a toy model consisting of a massive scalar field t) in a one-dimensional 
expanding background. In terms of the comoving coordinate x, the model is described by 
the Lagrangian 
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FIG. 1: (a) Energy density as a function of position for a typical final state: m = l,h = 0.5. (b) 
Snapshots of field profile <?!>(x,t) for a typical oscillon, taken at intervals of At = 0.1. 



where an overdot denotes a derivative with respect to t and a prime denotes a derivative 
with respect to x. Note that this potential has a unique minimum at = 0, with no other 
extrema or inflection points. It does not support static solitons, but does support oscillons. 
This Lagrangian leads to the equation of motion 

For simplicity we take the expansion rate H = d{t)/a{t) to be constant. Although we hope 
that our model can shed light on the reheating phase of inflation, during which H would 
be rapidly decreasing, for our purposes it is only important that the expansion rate be slow 
compared to the typical time scales of the oscillon. 

All known oscillon solutions have been found for massive fields, with frequency of oscil- 
lation below the threshold for fundamental excitations, so we have chosen a massive field. 
The sign of the 0^ term, the leading nonlinearity, is crucial: it has been chosen so that to 
leading order this term decreases the frequency of small oscillations. With the other sign we 
do not see any oscillons form. Finally, the 0^ term exists to eliminate instabilities at large 
field values introduced by the choice of sign of the 0^ term. Because the oscillons involve 
only moderate excitations of the field, this term does not play an important role in their 
dynamics. 

We have scaled the self-interaction terms by a small parameter g. The exact meaning 
of this parameter is somewhat subtle, since by defining a new field variable = ^/g(f) we 
can shift g completely outside the Lagrangian ([1]), which means that it no longer appears 
in the classical equations of motion. The full quantum theory is still sensitive to the value 
of g, but only in the combination %; therefore a classical approximation involving small 
h is equivalent to one involving small g, and the results of our classical analysis are valid 
provided that hg is small. ^ 



^ Also, throughout this paper the quantity we caU m has dimensions of [length] ^, not energy, and is thus 
really the mass of a single elementary quantum divided by h. 
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SMALL-AMPLITUDE ANALYSIS OF OSCILLONS 



To describe the oscillons we expect to emerge from the thermal background in our sim- 



ulation, we use a small-amplitude analysis, following |22j. In particular (e.g. see [l9|]) in a 
static background the equation of motion ([2]) supports solutions that are localized in space 
and periodic in time. These oscillons can be expanded in a small parameter < e -C 1. The 
amplitude of the oscillations scales with e, the spatial width of the oscillon is proportional to 
e~^, and the oscillation frequency is given hj u = m\/\ — . Numerical calculations indi- 
cate that there is an upper bound on how large e can be for a stable solution: < e < Cc < 1, 
where ~ 0.25. 

Here we extend these results to an inflationary background. A sufficiently fast expansion 
will result in a horizon that is comparable to the minimum oscillon width (proportional to 
e~^), preventing oscillons from ever forming. If the expansion rate is slow, long lived oscillons 
can occur in the range ^JlTpm ^ e ^ 1, as shown below. In this regime the oscillon radiates 
energy away (in the form of scalar field waves) at a rate that is exponentially small in the 
dimensionless ratio me^ /H. If this ratio becomes sufficiently large, numerical calculations 
show that the oscillon is not affected by the expansion in any way that we can detect. 



Separation of scales 



We work in static patch coordinates on de Sitter space (e.g., see [23i]) where the metric 
takes the form 

ds^ = -(1 - X^H^) dT^ + (1 - X^HY' dX^ . (3) 

These coordinates are valid for \XH\ < 1. In these coordinates the equation of motion ([2]) 
becomes 

_ ^2^2) <^ + 2^^'<^' - (1 - X'H')^" = -m' (0 - g<p' + g'^') . (4) 
Here an overdot indicates a derivative with respect to T and a prime indicates a deriva- 



tive with respect to X. Now we follow [l9j and change to variables x = mXe and 



r = mTvT^-^, where < e ^ 1 is the small parameter mentioned above. Note that 
since the horizon distance is and the oscillon width is 1/me, to obtain a stable oscillon 
solution we expect H -C me. We therefore let H = mHe^, where H is taken to be a small 
dimensionless number. The power of e in this expression is important: If instead we take H 
to be 0{e) then no oscillon solution is possible, while if we take H to be 0{e^) then (to the 
order we are working) the oscillon does not feel the expansion. 
In terms of x aiid r the equation of motion is 

+ 2x^^6^(0), - (1 - x'^e2)6'(0)xx = + 9'<f^') . (5) 

We expand (p in powers^ of e as 

0(X, ^) = bMx, r) + e^Mx, r) + e^Mx, r) + . ■ ] , (6) 



^ Note that the expansion ([6]) involves only odd powers of e. This is consistent with ([5]) because this equation 
is odd in 0, and depends on e via only. 
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and seek oscillons, solutions that are periodic (with period 27r) in r and locahzed in space. 
In fact, as we will show below, these solutions are neither strictly periodic nor localized, but 
these departures are exponentially small and occur over time scales much longer than the 
one given by r. 

We now substitute ([HD into the governing equation and solve it order by order in e. 
At 0(e), 

(0l)rr + 01 = O. (7) 

Thus 01 = /(x)^"*^ + where the profile f{x) remains to be determined. We can 

find this profile by considering the 0{e^) equation 

(03).. + 03 = (1 - X'i^')(0l)r. + (0l)xx + <PI ■ (») 

The key point here is that the only way that 03 can be periodic in r is if the forcing terms 
on the right-hand side of the equation are orthogonal to e^*"^, since otherwise 03 would have 
a component that grows linearly in r. Therefore we set the Fourier coefficients of e^*"^ to 0, 
and obtain a self-contained ordinary differential equation for the oscillon profile: 

+ {x'H' - 1) fix) + Vixfnx) = . (9) 

For a "perfect" oscillon, a localized (exponentially decaying as x ±cxd) spatial profile is 
needed. This, however, is not quite possible. The behavior of the above equation changes at 
X ~ ±l/-ff: for small |x| the equation is essentially the same as the equation for a flat-space 
oscillon, while for large |x| the x^H^ term dominates, giving oscillatory behavior that causes 
the oscillon to radiate (a small amount of) energy away. Note that this change in behavior 
occurs well before the horizon, which is at x = il/l^-^)- 



Asymptotic Behavior 

First we examine the regime where |x| -C l/-f^. Here the term arising from the expansion 
of the universe is negligible and the equation reduces to the case of a static background. 
Therefore, we can take 

1 /s" 1 

Note that since < if ^ 1, the oscillon behaves like ^J^e~^ in the region 1 <^ x <^ 

In particular, while its amplitude does not quite vanish (since we cannot take |x| ^ oo in 
the equation above), it fails to be fully localized only because of an exponentially small tail, 
as is shown below. 

We now relax the upper bound and consider the region where x ^ 1 {^"^ entirely similar 
calculation can be done where x ^ ■"!)• From the reasoning above, we see that on the left- 
hand side of this region, the field is exponentially small. Thus, provided that the field does 
not grow too much as x increases in this region, we are justified in neglecting the nonlinear 
term in the equation, though we must now take into account the effects of the expansion. 
Thus we obtain 

g + /(x^/^^-i)=o. (11) 
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This equation can be put into a more familiar form by defining a new coordinate y = xH. 
Then we obtain 



-H' 



,11 



y'f = -f. 



(12) 



This is exactly the Schrodinger equation for the wave function of a particle in an upside-down 
harmonic oscillator with energy E = —1, with the dimensionless Hubble constant H playing 
the same role as Planck's constant h in the analogous quantum mechanics problem. At small 
y the particle is in the "classically forbidden" region of the potential and the wavefunction 
is exponentially suppressed, corresponding to the exponential tail of the oscillon. For any 
nonzero value of H, we eventually enter a "classically allowed" region and the wavefunction 
becomes oscillatory, which for our classical field profile indicates an outgoing wave that car- 
ries energy away. Therefore this situation looks exactly like a quantum mechanics tunneling 
problem, and since we are interested in small H we can use a semiclassical approximation 
to solve it.^ 

Choosing an outgoing wave boundary condition as ?/ — oo, the standard WKB connection 
formulae 2J] give us the following relation between the wavefunctions on either end of the 
turning point at y = 1: 



A 



1 



,1/4 



exp 



H 



/2 



l)dy' + 



A 



2^1/4 



exp 



1 

'h 



^{l-y'^)dy' 



where the left-hand side is valid for y S> 1, the right-hand side is valid for ?/ -C 1, and A is 
an overall constant. 

Performing the integrals, keeping only the leading dependence, and replacing y with xH, 
we obtain 



/(x) ~ ^exp 
A 



I'm 



X 



exp 



for 1 < X < 1/H , and 
for x>l/-^- 



(13) 
(14) 



As expected, the small x behavior of fll3l) is of exactly the correct form to fit the large x 
asymptotic behavior of (ITO!) . Matching to this result sets A to exp (—2^) and hence fixes 

the coefficient of the outgoing wave ( fT4l) . 

Returning to our original variables and putting together the pieces, we find the following 
expressions for the oscillon 



4>{X, T) ^ eJ — cos (mT^/l — e^^ sech {mXe) for |X| <^ e/if , and 

V oQ V / 



0(X,T)^2e. 



m\x\g 



COS m 



iTVT 



1 



mXm 



(15) 

for e/H^\X\<l/H. 

(16) 



^ Of course, exact formal solutions for the wavefunction of an upside-down oscillator exist in terms of 
Hermite functions. However, extracting the asymptotic behavior from these special functions is somewhat 
messy; for our purposes we can obtain equivalent results simply using WKB. 
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From here it is easy to compute the relevant components of the stress tensor T^^ and find 
the rate of energy fiux. If we take the energy stored in a region R io he E = 
then using the conservation of the stress tensor and (fT6!) we obtain 



±E 

dT 



T\dR 



-(1 



32 
3^ 



exp 



Trme 



H 



(17) 



where the region R is taken to be a symmetric interval R = [—X, X] and X is far enough 
from the origin that (fT6!) holds. Note that to leading order all X dependence comes from 
the curvature of the metric; if we restrict attention to e/ H X <^ ^/H (i.e. we consider 
only a small neighborhood of the oscillon) then space looks almost fiat and we obtain 



dT 



32 2 3 

— m e exp 
3^ 



Time 
H 



(18) 



This is our main result. While an oscillon can live forever on a fiat background (at least as 
a formal perturbation series), this is no longer the case in a de Sitter universe; instead it is 
forced to radiate energy away, albeit through a mode that is exponentially suppressed. 



THERMAL INITIAL CONDITIONS 

We would like to start our simulations using initial conditions mimicking those of the 
interacting field theory defined by the Lagrangian ([T]) at nonzero temperature T. However, 
constructing this equilibrium is quite difficult; although we are essentially interested in 
classical physics, a classical treatment of this field theory at nonzero temperature suffers 
from the Jeans paradox. In a quantum treatment we avoid this problem, but we are still 
unable to systematically take into account the nonlinear terms in the Lagrangian: as is 
shown below, we will be interested in temperatures T ^ Tn/ which is precisely the regime 
where finite temperature perturbation theory fails. 

We shall thus take a different approach and generate our initial conditions to simulate 
thermal states of the free massive scalar field. We note that for the reasons stated above these 
quasi-thermal initial conditions are probably quite far from the true thermal equilibrium 
of the full interacting theory; hence the parameters T and h should be thought of more 
as measures of the amplitude and width of our distribution in momentum space than as 
directly physically relevant quantities. We show, however that provided T is sufficiently 
high, all other numerically feasible variations of these parameters produce oscillons in copious 
numbers, leading us to believe that our results are independent of any particular details of 
the initial conditions. 

To construct these conditions, we return to comoving coordinates. Since our real interest 
is in a numerical simulation we impose both infrared and ultraviolet cutoffs, placing the 
system in a box of comoving size L and on a regular lattice with spacing Ax. We replace 
the spatial derivatives by finite differences (see the Numerical Simulation section) and label 
the free field's normal modes by fc„ = 27rn/L, where n = —N/2 + 1 . . . N/2 and = L/Ax 
is the number of lattice points. Finally we take the scale factor at this time to be oq. 

On this lattice each free mode is described by a harmonic oscillator with frequency 

/2sin^ , , , 
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The initial conditions for are then given by 



N/2 

(a;^,t = 0)= ^ 

n=-N/2+l 
N/2 



h 



2LaQL0n 



[anc''^"'''' + a*e , and 



n=-N/2+l 



i V 2Lao 



(20) 



where a„ is a random complex variable with phase distributed uniformly on [0, 27r) and 
magnitude drawn from a Gaussian distribution such that 



huJr, 



- I coth ■ ^ 
2 V 2T 



(21) 



This is the usual amplitude distribution for a quantum harmonic oscillator [251] with the 
zero-point motion subtracted. On average, these initial conditions assign energy T to modes 
with hujn ^ T, in agreement with equipartition. The energy per mode goes rapidly to zero 
for fkVn ^ T, giving a total energy density that scales like T'^/h for T 3> mh, as usual for 
blackbody radiation in one dimension. Oscillons will form from the energy density in modes 
with wavelengths of order 1/m, which scales like mT. For oscillons to form, this energy 
density must be at least of order m^/g, so that the fields will have amplitude l/^/g and 
the nonlinear interaction terms can balance the dispersive gradient terms. Therefore we will 
need an initial temperature oiT ^ m/g to form oscillons. 

We have simply subtracted off the zero-point quantum fluctuations in the field. Although 
it is well-known that these fluctuations have significant consequences for the evolution of a 
classical field in an expanding background [26|], these effects are most important at very long 
length scales, of order 1/H. For small g, the oscillon solution comprises many fundamental 
excitations of the field at the length scales of order 1/m that are relevant to oscillon formation 
and stability. Thus the quantum effects we are neglecting should not change our results 
significantly. Alternatively, this quantum prescription can simply be thought of as describing 
a classical equilibrium with short- distance cutoff T/h. (A method to eliminate such cutoff 
dependence of classical simulations is discussed in |27|). 

As the universe expands, it cools and loses energy according to 



dE 
~dt 



— / p{x, t)<x{t)dx . 



(22) 



where the pressure density is 



+ 



/\2 



a{ty 



m 



2^ 3^ 



(23) 



In equilibrium at temperatures much greater than hm, if we neglect the interaction terms, 
the system looks like massless radiation, with pressure density approximately equal to its 
energy density. In equilibrium at low temperatures, the field is slowly varying in space, 
with only small values of k excited, and of small amplitude. In that case the gradient and 
nonlinear terms are negligible and 0^ ~ m^0, so the pressure goes to zero and the system 
behaves like pressureless dust. 
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NUMERICAL SIMULATION 



We discretize x at the level of the Lagrangian in Eq. ([T]), working in natural units where 
m = 1. For the space derivatives we use ordinary first-order differences, 

(24) 

where 0„ refers to the value of at lattice point n. We work on a regular lattice with 
spacing Ax and impose periodic boundary conditions. Varying this Lagrangian yields lattice 
equations of motion with second-order space derivatives, 

+ -^Mt) = a{t)^Ax^ ' 9Mt) - 9 Mt) ■ (25) 

We then express this second order equation as a set of coupled first order differential 
equations and step forward in time using a standard fourth-order Runge-Kutta integrator 
(e.g., see jl^). The Courant condition requires that we maintain At < a{t)Ax for numeric 
stability. Of course this allows us the possibility of rescaling At as the simulation runs; this 
rescaling stops when we reach a maximum value of At = 0.01. 

As the universe expands, the oscillons maintain a fixed size in physical units. On the 
other hand, our lattice expands with the universe. Thus we add new lattice points whenever 
the lattice spacing exceeds some fixed size Axmax in physical units. This is accomplished by 
refining the lattice, doubling the total number of lattice points and bringing the lattice spac- 
ing back to Axiaax/2 in physical units. We assign values to the field for the new intermediate 
lattice points by linear interpolation. 

We performed several checks on our numerics. We would like to take At sufficiently small 
that given a set of initial data any further reduction of At does not significantly change the 
final configuration after a run. In the strongly nonlinear regime in which we work, this turns 
out to be a technically very difficult goal to achieve, with different timesteps often resulting 
in significantly different final configurations. In fact, despite our best efforts the upper end 
of Fig. [21 where T > 6, suffers from this problem, so we cannot guarantee the validity of this 
region of the plot, although the remainder of the plots we present satisfy all of our tests. We 
also verify that all of our simulations maintain energy conservation (as given by Eq. fl22l) ) 
to better than one part in 10^, and that Axmax is small enough that further reduction does 
not significantly alter the final configuration obtained after a run. 

For Figure [2] we use an initial lattice spacing Axq = 0.0016 and an initial timestep of 
Ato = 0.0002. For Figure El we use Aa;o = 0.05/T and Ato = 0.0001, where T is the 
temperature of the run. For both we use Axmax = 1- 



RESULTS 

We simulate the scalar field in an expanding universe varying the initial temperature 
T and the value of h (which determines the effective coupling gh). We are interested in 
determining what fraction of energy in the universe ends up in oscillons; we estimate this 
as the integral of the energy density over regions of space at which the energy density is 
more than five times the average energy density, divided by the total energy. As shown in 
Fig. [H oscillons stand so much higher than the background fluctuations that this is a good 
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measure. For each run we expand until this quantity reaches a constant value, indicating 
that the system has stabilized. We then plot this quantity as a function of the parameters 
in the initial conditions in Fig. [2] and Fig. [31 

Fig. [T] shows a snapshot of the energy density as a function of position for a typical 
run. We see that the oscillons are very clearly defined, with energy densities far above the 
background from ordinary field fluctuations. At these late times, the ordinary fluctuations 
are of small amplitude, and have been redshifted to large spatial wavelengths; they have 
k ^ and uj ^ m. From Eq. (I23l) . we see that the pressure of such fluctuations vanishes: 
the 0^ and m^</)^ terms cancel, and the gradient and nonlinear terms can be neglected 
because the ordinary fluctuations are slowly varying in space and have small amplitude. So 
the total energy in ordinary fluctuations remains constant, corresponding to a pressureless 
dust of ordinary particles of mass km at rest. As the universe continues to expand, the 
energy density in these fluctuations decreases, since it scales inversely with the volume of 
the universe to keep the total constant. The oscillons, meanwhile, maintain fixed physical 
size, energy density, and total energy — and therefore zero pressure — but are brought to 
rest relative to the expanding background. As a result, they also appear as a pressureless 
dust at late times. 

In Fig. [2] we vary h in the initial conditions and examine the variation of the fraction 
of energy in the universe in oscillons. Each point on the graph is an average over many 
runs, with error bars indicating the standard deviation. (Because our initial conditions are 
random individual representatives drawn from the initial thermal distribution, our results 
can vary from one run to the next.) As we move to the right on the logarithmic scale on 
the horizontal axis, h becomes smaller and the quantum effects we have neglected are less 
important. We see that the fraction of energy in oscillons remains substantial, decreasing 
only gradually as h decreases to small values. 

Fig. E] shows the fraction of energy in oscillons for a range of initial temperatures. For T 
large enough, we see that the result saturates, so that higher initial temperatures no longer 
affect the final result. In the saturated regime, the system is just undergoing ordinary 
cooling, remaining at equilibrium as the temperature begins to decrease. It is only once 
the system cools below T ^ fn/g that the oscillons begin to emerge. So while we are 
always imagining our system starts out at very high temperatures compared to the energy 
scales relevant to oscillon formation, in practice we only need to start our simulations at 
temperatures just above the saturation point. 

DISCUSSION AND CONCLUSIONS 

We have considered a scalar field theory that supports stable oscillons. We place this 
theory in an expanding background and show that the oscillons are no longer completely 
stable but instead lose energy at a rate that is exponentially small in the size of the horizon. 
We also find numerically that quasi-thermal initial conditions result in the eventual forma- 
tion of oscillons in large numbers. Though the oscillons are large, coherent objects, they 
nevertheless form easily from a random superposition of momentum modes, suggesting that 
in some sense they are attractors in the solution space of the equations of motion. In our 
model they capture a significant portion of the energy of the universe, on the order of 50%. 
Though this paper deals with a one-dimensional case, we have seen qualitatively similar 
results in two dimensions. If this phenomenon persists in realistic models (and in a realistic 
number of dimensions), oscillons may have cosmological consequences, as discussed in Ref. 
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FIG. 2: Fraction of energy in oscillons, as a function of log(^). Simulation parameters are T = 2, 
g = 1, initial universe size Lq = 20, expansion factor = 2000, Hubble constant H = 0.02. 
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FIG. 3: Fraction of energy in oscillons, as a function of initial temperature. Simulation parameters 
are h = 0.5, g = 1, initial universe size Lq = 40/r, expansion factor = 40000/Loi Hubble constant 
H = 0.02. 
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